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SYSTEM AND METHOD FOR DETERMINING A DIELECTRIC PROPERTY 
ASSOCIATED WITH A SUBSTRATE 

CROSS-REFERENCE TO RELATED APPLICATION (S) 

[0001] This application discloses subject matter related 
to the subject matter disclosed in the following commonly 
owned co-pending patent application (s) : (i) ''SYSTEM AND 
METHOD FOR SPECIFYING A DIELECTRIC PROPERTY OF A SUBSTRATE 

UNDER DESIGN," filed , Application No.: 

(Docket Number 200312 073-1) , in the names 

of: Karl Joseph Bois, David W. Quint and Timothy L. Michalka; 
and (ii) ^^SYSTEM AND METHOD FOR DETERMINING A DIELECTRIC 
PROPERTY ASSOCIATED WITH A CONSTITUENT MATERIAL OF A 

COMPOSITE SUBSTRATE," filed , Application 

No.: (Docket Number 200312075-1), in the 

names of: Karl Joseph Bois, David W. Quint, Timothy L. 
Michalka and Peter Shaw Moldauer, both of which are hereby- 
incorporated by reference for all purposes. 

BACKGROUND 

[0002] In recent years, the computer industry has strived 

to introduce system links that operate in the gigahertz 
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regime . Due to the increase in data bandwidth and decrease 
in voltage margin, the attenuation and distortion caused by 
the channel is of much more concern than in past system 
generations. Hence, for successful modeling and design of 
signal interconnects, it is critical to determine the 
frequency variations of the dielectric characteristics of the 
signal channel, especially their frequency-dependent losses. 
As such, the accurate knowledge of electric properties of 
chip packaging and printed circuit board (PCB) materials is 
required - 

[0003] For example, in PCBs, a widely used material is 
what is commonly referred to as Fire Retardant (FR) -4 
material, which is relatively less expensive. FR4 is a 
material formed by glass strands embedded in an epoxy resin 
binder. Counter to its low cost benefit, the material 
exhibits noticeable attenuation at higher frequencies, for 
example at frequencies greater than 1 GHz. These losses are 
associated with the loss tangent of the material, wherein the 
complex dielectric constant, e^, of the material varies by 
frequency and may be expressed by the following equation: 

= ^'r - 3^"r/ wherein 

e'r is the relative permittivity of the material 
which varies as a function of frequency; and 

is the loss factor of the material which varies 
as a function of frequency. 
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Loss tangent, tan 5, may be defined by the following equation 
in terms of the relative permittivity, e'r/ B-nd the loss 
factor, 

tan 5 = £"r/£'r 

[0004] Unfortunately, when modeling transmission lines 
embedded in FR4, data supplied by most manufacturers for this 
parameter is usually measured at relatively low frequencies, 
for example, 60 Hz or 100 MHz. On the other hand, existing 
measurement techniques such as coaxial or waveguide 
techniques are not only impractical for In-sltu measurements, 
they are generally inadequate for capturing higher- order 
effects. Hence, a method for measuring the dielectric 
properties, such as the dielectric constant and loss tangent, 
of this material and other materials in the gigahertz regime 
is needed - 

SUMMARY 

[0005] A system and method are disclosed that provide for 
determining the dielectric properties associated with a 
substrate . In one embodiment , a network analyzer measures 
scattering parameters for at least two lines of substantially 
identical cross-section embedded within the substrate over 
a specified frequency range. A first engine determines a 
complex propagation constant based on the scattering 
parameters and defines the complex propagation constant in 
terms of an attenuation component and a phase component. A 
second engine, responsive to the phase component, determines 
a relative permittivity parameter associated with the 
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substrate over the specified frequency range. A third 
engine, responsive to the attenuation component and the 
relative permittivity parameter, performs a least squares 
analysis to determine a loss tangent parameter associated 
with the substrate over the specified frequency range. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0006] FIG. 1 depicts a schematic diagram of one 
embodiment of a system for analyzing a substrate; 
[0007] FIG. 2 depicts a schematic diagram of a pair of 
probes for measuring scattering parameters of a stripline via 
plated through holes (PTHs) ; 

[0008] FIG. 3 depicts a functional block diagram of one 
embodiment of a system for determining the dielectric 
properties associated with a substrate; 

[0009] FIG. 4 depicts a flow chart of one embodiment of 

a method for determining the dielectric properties associated 
with a substrate; 

[0010] FIG. 5 depicts a graph of attenuation component as 

a function of frequency; 

[0011] FIG. 6 depicts a graph of relative permittivity and 

loss tangent as a function of frequency; 

[0012] FIG. 7 depicts a graph of relative permittivity as 
a function of frequency for several types of FR4 -based 
printed circuit boards (PCBs) ; and 

[0013] FIG. 8 depicts a graph of loss tangent as a 
function of frequency for several types of FR4 PCBs. 
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DETAILED DESCRIPTION OF THE DRAWINGS 

[0014] In the drawings, like or similar elements are 
designated with identical reference numerals throughout the 
several views thereof, and the various elements depicted are 
not necessarily drawn to scale. Referring now to FIG. 1, 
therein is depicted an embodiment of a system 10 0 for 
analyzing a substrate, which is illustrated as a printed 
circuit board (PCB) 102. A work surface 104 supports the PCB 
102 and provides a physical area where an engineer can 
interact with the PCB 102. A microscope assembly 106 is 
utilized by the engineer to enlarge the minute elements, such 
as plated through holes (PTHs) and microstrips, of the PCB 
102 by way of one or more optical lenses. A network analyzer 
108 derives the electrical properties of the PCB 102 using 
probes 110 and 112 which provide coupling to the PCB 102. 
In particular, probes 110 and 112 inject or extract energy 
relative to the PCB 10.2 in order to measure various 
electrical parameters such as a standing wave-ratio or 
scattering parameters. A computer 114 is associated with 
network analyzer 108, or may be integrated therewith, to 
receive and process the measurements and data collected by 
network analyzer 108. A monitor 116 is positioned in 
communication with the computer 114 and presents the raw and 
processed measurements to the engineer via a display. 
[0015] In operation, as will be explained in more detail 
hereinbelow, the probes 110, 112 coupled to the network 
analyzer 108 enable in-situ scattering parameter measurements 
for at least two lines of substantially identical cross- 
section embedded within the PCB 102 over a specified 
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frequency range. The network analyzer 108, responsive to the 
in- situ measured scattering parameters, determines the 
relative permittivity, and loss tangent, tan 5, 

associated with the PCB 102 by utilizing a complex 
propagation constant and a cascade representation of uniform 
transmission to model the PCB 102. The loss tangent, in 
turn, may be employed to determine the complex permittivity, 
e^/ of the material, for example. 

[0016] FIG- 2 depicts a system 2 00 for the measuring of 

scattering parameters in a substrate which is illustrated as 
PCB 202, which includes strata 204, 206, 208, and 210 formed 
of a fire-resistant cured heterogeneous material which may 
comprise a composition of glass strands and resin binder. 
In one embodiment, the PCB 202 may be a fire retardant (FR) -4 
material of any grade, such as Nelco 4000-6, 4000-13, or 
4000-13 SI, for example. Stripline 212 is positioned in the 
layer 208 and PTHs 214 and 216 provide surface contacts for 
the stripline 212. Similarly, stripline 218 is positioned 
in the layer 2 08 and PTHs 22 0 and 222 provide surface 
contacts for the stripline 218. It should be appreciated 
that although PTHs are depicted in FIG. 2, other types of 
connectors may be employed such as subminiature version A 

(SMA) connectors. As illustrated, the striplines 212 and 218 
are positioned within the same layer of the PCB 202 and have 
substantially identical cross-sections, but have different 
lengths, e.g., stripline 218 is longer than stripline 212. 
It should be further appreciated that although the line is 
described as a stripline, other types of transmission lines, 
including non-embedded conductors such as microstrips, for 
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example, may be selected for the determination of dielectric 
properties of the substrate • 

[0017] Once the network analyzer is calibrated, probes 224 
and 226 are coupled to PTHs 214 and 216, respectively, in 
order to perform in-situ scattering measurements of the PCB 
202. For the purposes of measurements and calculations, the 
stripline 212 and its PTHs 214 and 216 are considered as a 
uniform transmission line which is attached to connectors 
with unknown or not easily extractable electrical 
characteristics. This arrangement is modeled as a cascade 
representation of a uniform transmission line between two 
connectors for the length, 1, of stripline 212 (see, for 
example, Janezic, M.D., and J. A. Jargon, ^^Complex 
Permittivity Determination from Propagation Constant 
Measurements," IEEE Microwave and Guided Letters, vol. 9, no. 
2 (1999), pp. 76-78) , which may be described by the following 
equation : 

= XT^Y, wherein 

is the scattering cascade matrix, i.e., a square 
array of numbers consisting of the transmission 
and reflection coefficients associated with the 
two connectors, i.e., PTHs 214 and 216, and the 
transmission line, i.e., stripline 212; 
X and Y represent the connections PTH 214 and 216, 
respectively, from the calibrated plane of the 
network analyzer to the stripline 212 ; and 
T^ represents the homogenous transmission line, 
i.e., stripline 212. 
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[0018] As will be described in more detail hereinbelow, 
by collecting in-situ scattering parameter measurements of 
the PCB, utilizing a cascade representation of uniform 
transmission through the conductor line, e.g., the stripline, 
a complex propagation constant may be computed for the line 
conductor, which will be employed by the network analyzer to 
determine the dielectric properties associated with the 
substrate. Accordingly, it should be appreciated that the 
scheme described herein enables the determination of the 
dielectric properties of post-assembly substrates without 
destroying the substrate. Moreover, this scheme employs 
embedded transmission lines and does not necessarily rely on 
any explicit modeling of launch points, e.g., PTHs . Further, 
it will be seen below that the teachings described herein 
permit the differentiation between conductor and dielectric 
losses, after taking into account such higher-order effects 
as the skin effect. 

[0019] FIG. 3 depicts a functional block diagram of one 
embodiment of a system 300 for determining the dielectric 
properties associated with a substrate. Scattering parameter 

(S-parameter) measurements 302 for at least two lines of 
substantially identical cross-section embedded within the 
substrate over a specified frequency range are provided to 
a complex propagation constant engine 304, i.e., a first 
engine. In one embodiment, the frequency range is from about 
300 KHz to about 8.5 GHz. The complex propagation constant 
engine 3 04 determines a complex propagation constant based 
on the scattering parameters using the cascade representation 
of a uniform transmission line described hereinabove. In 
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particular, as will be evident to those skilled in the art, 
the complex propagation constant of a uniform transmission 
line is determined by ^'de- embedding" the connection points 
that are attached to it. The complex propagation constant 
engine 304 then defines the complex propagation constant in 
terms of an attenuation component 306 and a phase component 
308, both of which are frequency dependent, i.e., in terms 
of a real part and an imaginary part. More specifically, the 
complex propagation constant, y, which varies as a function 
of frequency, is defined as follows: 

Y(f) = Of(f) + j3(f), wherein 

a(f) is the frequency dependent attenuation 
component (i.e., the real part); and 
j3(f) is the phase component, which is the 
frequency dependent imaginary portion of yif) . 

[0020] A relative permittivity constant engine 310, i.e., 

a second engine, responsive to the phase component 3 08, 
determines the relative permittivity 312 associated with the 
substrate over the specified frequency range. More 
specifically, the phase component, P(f), may defined as 
follows : 
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3(f) = [CO /c, wherein 

CO is the angular frequency; 

c is the speed of light, i.e., approximately 
2.998 X 10®m/s; and 

s'r is the relative permittivity of the substrate. 

By defining the complex propagation constant in terms of the 
attenuation component 3 06 and the phase component 3 08, the 
relative permittivity may be determined as a function of the 
frequency since the speed of light and the phase component 
are known . 

[0021] A least squares engine 314, i.e., a third engine, 
responsive to the attenuation component 306 and the relative 
permittivity 312, performs a least squares analysis to 
determine the loss tangent associated with the substrate over 
the specified frequency range. More specifically, the 
attenuation component, a, may be decomposed into attenuation 
attributable to the conductor and attenuation attributable 
to the substrate as illustrated by the following ecjuation: 

a = + Ofd/ wherein 

is the attenuation attributable to the 

conductor, i.e., the stripline or microstrip; and 

is the attenuation attributable to the 
dielectric material . 



[0022] The attenuation attributable to the conductor, a^, 
may be further defined in terms of the line resistance and 
the skin effect, which is the tendency of alternating 
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currents to flow near the surface of a conductor thus being 
restricted to a small part of the total sectional area and 
producing the effect of increasing the resistance, as 
follows : 

= [Rdc + Rsi/'f]/2Zo, wherein 
Rdc is the line resistance; 

Rg is the resistance due to the skin effect; and 
Zq is the characteristic impedance. 

[0023] By utilizing this model, the least squares engines 
314 may determine the attenuation attributable to the 
conductor, oc^, since the line resistance, R^c/ the resistance 
due to the skin effect, Rg, and the characteristic impedance, 
Zq, are constants. The attenuation attributable to the 
attributable to the dielectric material, a^, may further be 
defined in terms of relative permittivity and the loss 
tangent as follows: 

= [nf/c] (tan 5)-/^'^ 

[0024] As the line resistance, R^c/ and the resistance due 

to the skin effect, Rg, can be treated as constant parameters 
and, at low frequencies, the loss tangent, tan 5, is nearly 
constant, the attenuation component, a, may be written in 
terms of the attenuation attributable to the conductor, a^, 
and the attenuation attributable to the dielectric material, 
a^, as a polynomial function such that: 
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a = [l/2Zo]Rdc + [l/2Zo]Rs"/f + [n/c] (tan 6)^^E\'f, 

wherein [l/2Zo]Rdc; [l/2Zo]Rs; and [n/c] (tan 
are constants. 

From this equation, a least squares model may be set up as 
follows : 

X(f) = [1 /f f] 
b = (X'X) X'a 

b = [Rcic/2Zo Rs/2Zo constant] 

Since is known, all other quantities may be determined. 
In particular, the loss tangent, tan 5, may be determined, 
which in turn can be employed in further calculations such 
as determining the substrate's complex dielectric constant . 
[0025] FIG- 4 depicts a flow chart of one embodiment of 
a method for determining the dielectric properties such as 
the loss tangent associated with a substrate. At blocJc 400, 
scattering parameters are measured for at least two lines of 
substantially identical cross-section embedded within the 
substrate over a specified frequency range. In one 
implementation, the range is from about 300 KHz to about 8.5 
GHz. Additionally, in a further embodiment, the methodology 
measures scattering parameters for two transmission lines 
that have similar launch points, e.g., both have SMA 
connectors, and different lengths. At bloclc 402, a complex 
propagation constant is defined based on the scattering 
parameters. At block 404, the complex propagation constant 
is defined in terms of an attenuation component and a phase 
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component. At block 406, based on the phase component, the 
relative permittivity associated with the substrate over the 
specified frequency range is determined. At block 4 08, based 
on the attenuation component and the relative permittivity, 
a least squares analysis is performed to determine the loss 
tangent associated with the substrate over the specified 
frequency range wherein the conductor and dielectric losses 
are differentiated. Accordingly, this methodology may be 
employed in conjunction with any test vehicle for monitoring 
material properties of composite substrates. For example, 
the instant methodology presented herein may be employed for 
PCBs, chip packages, or silicon, for example. Furthermore, 
the various engines that may be utilized for determining the 
dielectric properties of the material can be implemented in 
software, hardware, firmware, or any combination thereof 
supported by a computer environment such as the system shown 
in FIG. 1. 

[0026] In an experimental verification of the teachings 
presented herein, three boards were produced with the 
following specifications : 

o five striplines of 2", 4", 12", and 20" in length; 
o three types of FR4 : Nelco 4000-6, 4000-13, and 
4000-13 SI; 

o all striplines routed on the same layer; and 
o SMA connectors . 

The cross-sectional specifications for stripline width, w, 
substrate dielectric thickness, h, and stripline thickness. 
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t, for the three boards are presented in the following table. 
Table 1. 

Table 1. Cross -Sectional Specifications 



Material 


w 


h 


t 


4000-6 


584 um 
(23 mils) 


660 um 
(26 mils) 


3 5.6 jam 
(1.4 mils) 


4000-13 


584 um 
(23 mils) 


533 um 
(21 mils) 


3 5.6 \im 
(1.4 mils) 


4000-13 SI 


584 lam 
(23 mils) 


521 lam 
(20.5 mils) 


3 5.6 \im 
(1.4 mils) 



[0027] FIG. 5 depicts a graph 500 of attenuation 
component, oc, as a function of frequency for the PCB produced 
with the aforementioned Nelco 4000-6 material. As expected, 
the contribution of dielectric loss to the attenuation 
component is linear with frequency and that of conductor loss 
is proportional to the square root of frequency, f. 
Additionally, the agreement between the measured and 
reconstructed values of the attenuation component, a, is very 
good. This indicates that the assumptions used for fitting 
the data according to the least squares model are indeed 
valid. 

[0028] FIG- 6 depicts a graph 600 of the relative 
permittivity, s'r/ loss tangent, tan 5, as a function of 

frequency for the PCB constructed with the 4000-6 material. 
In particular, the graph 600 illustrates cases for when the 
contribution of conductor losses are removed from the total 
losses and when they are not. As illustrated in graph 600, 
when the skin effect losses are not subtracted from the total 
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losses, the measurement of the loss tangent, tan 5, is much 
higher and also presents frequency characteristics, which 
could only follow a Kramers -Kronig relationship if the 
material under test, i.e., the Nelco 4000-6 material, was a 
conductor (which is not the case) . Additionally, graph 600 
illustrates that the loss tangent, tan 5, is fairly constant 
over a wide range of frequencies, which further validates the 
systems and methods described herein. As a further 
demonstration of the effectiveness of the teachings presented 
herein, it is interesting to note from FIG. 6 that the value 
of loss tangent for the PCB constructed with the Nelco 4 000-6 
is around 0.02 at 1 GHz, which is the quoted value provided 
from the manufacturer of the Nelco 4000-6 material. 
[0029] FIG. 7 depicts a graph 700 of relative permittivity 
as a function of frequency for the three types of FR4 PCBs . 
Similarly, FIG. 8 depicts a graph 800 of loss tangent as a 
function of frequency for the same three boards. As 
illustrated, with reference to FIGS. 7 and 8, for frequencies 
below 1 GHz, the variation of e'^. a function of frequency 
exhibits behaviors that are similar to the observed tan 5, 
when the effect of conductor losses was not included as part 
of the loss mechanism. In the frequencies of interest for 
high speed system links, i.e., greater than 1 GHz, the error 
is minimal as the contribution of internal inductance is very 
small. In this region, the variation of E \ as a function of 
frequency is almost constant. 

[0030] Accordingly, a simple but accurate method for 
determining the complex permittivity and associated values 
of PCB materials has been presented that employs a specific 
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formulation for determining the complex propagation constant 
of a uniform transmission line by de-embedding the connection 
points that are attached to it. By modeling the physical 
behaviors of dielectric materials and conductors, the 
conductor and dielectric losses are separated from each 
other. This permits the extraction of relative permittivity 
and loss tangent of a variety of substrates , including 
composite substrate materials. 

[0031] Although the invention has been particularly 
described with reference to certain illustrations, it is to 
be understood that the forms of the invention shown and 
described are to be treated as exemplary embodiments only. 
Various changes, substitutions and modifications can be 
realized without departing from the spirit and scope of the 
invention as defined by the appended claims. 
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